Hilbert fluid dynamics equations expressed in 
Chapman-Enskog pressure tensor and heat 

current 

Soderholm, Lars H. 

Mekanik, KTH, SE-10044 Stockholm, Sweden 
lars . soderholmOmech . kth . se 



1 Abstract 

The connection between the Chapman-Enskog and Hilbert expansions is investigated 
in detail. In particular the fluid dynamics equations of any order in the Hilbert 
expansion are given in terms of the pressure tensor and heat current of the Chapman- 
Enskog expansion 
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2 Introduction 

The most convenient way to calculate heat current and pressure tensor is that of the 
Chapman-Enskog method, see Grad [4]. In particular, the terms up to second and 
higher order have been calculated and at least to second order there is an agreed 
convention of notation for the terms. But in many occasions the Hilbert method 
is preferred because of the difficulties with the Burnett equations, see Sone [6]. 
Already in the Handbuch article [4], Grad made clear that the Chapman-Enskog 
method and Hilbert method are two ways of expressing the same thing. But Grad 
adds "to confirm this directly is an intricate task due to the different formalisms". 
So, his proof is somewhat abstract. It is the object of the present work to do the 
explicit confirmation that Grad does not do. 

The benefit of this is that the fluid dynamics equations of the Hilbert expansion 
are given explicitly in terms of the pressure tensor and heat current of the Chapman- 
Enskog method. This means that results that have been derived by the Chapman- 
Enskog method directly can be used in the Hilbert expansion. In the interesting 
paper [3] by Chekmarev and Chekmareva multiple scale methods are applied in the 



Hilbert expansion. The results of the present paper should considerably simplify and 
generalize their results to second order. In the paper [5] the results by Chekmarev 
and Chekmareva were generalized but derived from the Burnett equations. This 
could be considered unsatisfactory as the Burnett equations are mathematically 
unstable. But here we have been able to show that this result was nevertheless 
correct. We have shown that the fluid dynamics equations of the Hilbert expansion 
in fact can be obtained from the corresponding equations of the Chapman-Enskog 
expansion by expanding the fluid dynamics fields. 



3 Chapman-Enskog expansion and expansion of fluid 
dynamics variables 

We write the Boltzmann equation in dimensionless variables (e is the Knudsen num- 
ber) 

eVf = Q(fJ), 
where the streaming operator is (F is the force per mass) 

^ 9 9 T. 9 
V = h c h F . 

dt <9x dc 
See Cercignani [TjThe fluid dynamic fields are given by 

771 / = P ' 

m J fcd 3 c = pv=j, (3.1) 
m f 2 3 v 2 3k B T 

-Jj fcdc = " ( T + ^T ) = e 

In the Chapman-Enskog expansion is based on the important auxiliary condition, 
that only the zero order distribution function contributes to the flud dynamics fields. 
Hence (|3. 1 j) is satisfied also with / replaced by /oEoWe summarize the fluid dy- 
namical fields 

P 



w 



v 

T 



The components of w are written w a where a = 1...5. The Chapman-Enskog series 
is written (e is the Knudsen number) 

Ice = fcEo(w,c) +ef C Ei(w,c) + ... 

Here, fcEr(w,c) stands for a function of w and its r first derivatives and of c, the 
molecular velocity. In the sequel we don't write out c explicitly. 

Let us now assume that the fluid dynamics variables can be expanded in the 
Knudsen number 

w = wq + e^i + ... 
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Applying this expansion in fcEr{w,c) we find 

fcEr(w) = fcEr,o( w o) + efcEr,l( w 0, w l) + - 

Here fcEr,s is a function of wo,w\, ...w s and their r first derivatives. 

We first of all consider the zero order term fcEO- When we expand it we obtain 
to lowest order fcEO,o(wo). This is a Maxwellian for the fluid dynamics fields wq. 
In other words, this is the zero order term in the Hilbert expansion 

fcEo,o = fm- 

Let me write the sum of the first r terms in the Chapman-Enskog expansion and 
introduce this expansion of the fluid dynamics fields 

fcEO + tfcEl + ■■■£ n fcEr 

r 

= fcEo,o + e(fcEo,i + fcEi,o) + •••e r fcEr-s,s + 0(e r+1 ) 

s=0 

Thus we obtain 

p = m J fcEod 3 c = m J (fcEO + efcEi + • ••£ r fcEr)d 3 c 
= m J fcEO,od 3 c + em J (fcE0,i + fcEi,o)d 3 c + ... 

+e r m f j2fcEr-s,sd 3 c + 0(e s+1 ) 

Identifying the terms we find 



p r = m / ^ fcEr-s,sd 3 c. (3.2) 
In the same way we obtain 

jr = {py) r = ^2 PrsVs = m fcEr-s, S cd 3 C. (3.3) 

Finally we have 



s=0 J s=0 



v 2 3k B T 



Pr-s{2^2^s-p^p + —T s ) 

y [iTfcEr-s,sC 2 d 3 c (3.4) 



s=0 
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From (j3.2H3.4p we conclude that if we assume that the fluid dynamics field are 
expanded in a series in the Knudsen number, the corresponding terms to order r of 
the fields are given as the projections of 



y^JcEr-s,s- (3.5) 



s=0 

r 

We shall see that ^ fcEr-ss is m fact the coefficient /#> of e r in the Hilbert 

s=0 

expansion. 

4 The Hilbert distribution function expressed in the 
Chapman-Enskog distribution function 

Let us now consider the terms up to order r in the Chapman-Enskog expression of 
the collision term. We then expand the fluid dynamics fields and obtain 

r 

£ e p+q Q(fc Ep JcE q ) 

p+q=0 

r 

£ e p+q+i+j Q(fc Ep ,i, fcEtf) + 0(e r+1 ) 

p+q+i+j=0 

In the coefficient of e r we put p + i = k, q + j = I so that it becomes 

£ Q{fcEk-p,p, fcEl-q,q) 
k+l=r 

k I 

= £ QC/j fCEk-p,p, / j fCEl-q,q) 
k+l=r p=0 9=0 



We obtains sums of the kind 

If we do the corresponding thing in the streaming term we find that the coefficient 
of e r_1 is 

r-l 

£ fcEr~l-p,p 
p=0 

So from the Chapman-Enskog equations of order 0, 1, ...r we find that 

r r 

Q{^2 fcEr-p,p, fCEOfi) + Q(fCE0,0, £ fcEr-p,p) 
p=0 p=0 
r— 1 r— 1 

= - £ Q(%2 fCEk-p,p,^2 fCEl-q,q) (4.6) 
fc+Z=r p=l q=l 
r-l 

fCEr-l-p,p 

p=0 
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Now we proceed by induction. We assume that the sum (|3,5p to order r — 1 is 
the corresponding term in the Hilbert expansion 



r-1 



}ht-1 = y] fcEr-l-s,s- (4.7) 

We can then write (14.61) as 



QC/; fcEr-p,p, fm) + Q(fm, fcEr-p,p) 
p=0 p=0 

r-1 

= - QifHr-s, fHs) + DfHr-1 



s=l 



This means that 



satisfies the equation 



ic£r- s , s (4.8) 



s=0 



Q(f Hr jHo) + Q{fmJ H r) (4.9) 



r-1 



= - Qifnr-s, Jhs) + T^fm-i 

8=1 

for in the Hilbert expansion. We assume the existence of a solution and its 
uniqueness up to a linear combination of l,c,c 2 . But according to (|3.2H3.4|) the 
function (|4.8p also has the correct projections for /#> on the functions l,c,c 2 . So by 
induction this is the function fg r of the Hilbert expansion. We conclude that the 
simple relation between the Chapman-Enskog and Hilbert distribution functions is 

r 

fHr = J2fcEr-s,8- (4.10) 



5 Conservation laws and the Chapman-Enskog equa- 
tions 

Let us to start by writing down the conservation laws. They are 

P,t + V-j = 

j,t + V-J = pF (5.11) 
e,f + V-j e = F-j 

Here 

J =m j fccd 3 c = P+pvv = P + jv, (5.12) 
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where 

P =m 



J fc'c'd 3 c (5.13) 



and 

„2 



j e = f^fcd'c (5.14) 
= ev + q + P • v, 

q=^| c' 2 fc'd 3 c. (5.15) 
If we instead calculate the moments from fcEr we obtain to zero order 

3cEQ=m j fcEoccd 3 c = Pc£0+jv, (5.16) 

P C E0=m J f C Eoc'c'd 3 c = pl (5.17) 



ieCEO 



For r > 1 we have 



= J ^fcEocd^c (5.18) 
= ev + q C Bo+ p C£;o-v, 

qcso=f | c /2 /c'd 3 c = 0. (5.19) 

JcEr=™ y fcErCcd 3 C = PcEr, (5.20) 

P C Er=m J f C ErC'c'd 3 C (5.21) 

-^-fcErCd 3 c (5.22) 

= q C Er+~PcEr-V, 

<lCEr = - j c' 2 f C ErC'd 3 C. (5.23) 

Peer an d qcfir an d hence JcEr and j e CEr are functions of the fluid dynamics fields 
w and its r first derivatives. This means that we can expand the fluid dynamics 
fields in them to obtain 

PcEr = Pc£r,0 + ePc£r,l + ■■; 

where PcErs is a function of wq,wi, ...w s and its r first derivatives. 
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The Chapman-Enskog equations to order r are obtained if we make the approx- 
imations 



P -> ^£ s P CEs , (5.24) 

s=0 
r 

s=0 

in the expressions for J and j e and substitute into the equations of balance. The 
equations are then 

P.t + V-j = 

r 

j, t + V-(^e s J C £s) = pF (5.25) 

s=0 
r 



e ,t + v - (^e'jecss) = F j 

s=0 

6 Integrability conditions for the Hilbert expansion 

We consider the equation (|4.9p of the Hilbert expansion where we now replace r by 
r + 1. The condition of integrability is that is 

m \ 

mc \ Vf Hr d 3 c = 0. (6.26) 
mc 2 /2 J 

We now use the relation (|4,1U|) to express this condition in terms of the Chapman- 
Enskog distribution function. To that end we need the corresponding integrals for 
the Chapman-Enskog terms. To zero order we obtain 



m j Vf C Eod 3 c = /Ot + V-j, 
m J Vfc Er cd 3 c = it + V ■ Jceo-pF, (6.27) 



m / -J V fcEocd 3 c = e, t + V • j eCE o - F • j. 
To order r > 1 there will be no contributions to the time derivatives 

m [ Vf CEr d 3 c = 0, 



m J Vf CE rCd 3 c = V • Jcet, (6.28) 



,-,2 



m / —VfcErCd^C = V -ieCEv 
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If we now expand the fluid dynamics fields we obtain 
m J T>f C E0,sd 3 c = p Sjt + V-j s , 
m J Vf C E0, s cd 3 c = j a>t + V • J C eo,s-PsF, (6.29) 
m / ^fcEo^c = e s , t + V • j eCEOtS - F • j s . 
For r > 1 we have 

m J T>f C Er, s d 3 c = 0, 
m J VfcEr,sCd 3 c = V • J CEr ,s, (6.30) 



m 



2 



VfcErcd 3 c = V-j eCEryS . 



Now, using the expression for /#> in terms of fcEr-s,s we finally obtain the condi- 
tions of integrability to order r 

Pr,t + V-j r = 0, 

r 

Jr,t + V- ^ J C i?r- S , S -prF= 0, (6.31) 

s=0 

r 

e r ,t + V- jeOEOr-s.s ~ F • j r = 0. 
s=0 

They are the fluid dynamics equations of the Hilbert expansion. Note that the same 
equations are obtained if we start with the Chapman-Enskog equations of order r, 
expand the fluid dynamics fields and pick out the coefficient of e r . 

Let us now as an example consider the linearized Boltzmann equation and find 
the Hilbert fluid dynamics equations to second order. We have, see Chapman & 
Cowling 

Pcei = -2/x < Vv >, 

2 2 

P C E2 = ^(w 3 -w 2 )(VVT)-w 2 ^(VVp), (6.32) 
qcei = -kVT, 

qc,2 = 7[(y)Av+(^-^)V(V.v)]. 
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Hence, the Hilbert fluid dynamics equations to second order are 

P2,t + PoV ■ v 2 = 0, 
Pov 2 ,t = -V P2 + 2/x (Avi + |v(V • Vl )) (6.33) 

-[^(^3 - w 2 )] AVT + [w2t^]oAV Po , (6.34) 
^p( r 2,t + T V • v 2 ) = kqATi - [|^(0 4 - 9 2 )} A(V ■ v ). 
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